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Abstract. Holonomy corrections to scalar perturbations are investigated in the 
loop quantum cosmology framework. Due to the effective approach, modifications of 
the algebra of constraints generically lead to anomalies. In order to remove those 
anomalies, counter-terms are introduced. We find a way to explicitly fulfill the 
conditions for anomaly freedom and we give explicit expressions for the counter-terms. 
Surprisingly, the /i-scheme naturally arises in this procedure. The gauge invariant 
variables arc found and equations of motion for the anomaly-free scalar perturbations 
are derived. Finally, some cosmological consequences are discussed qualitatively. 
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1. Introduction 

Loop Quantum Gravity (LQG) is a tentative non-perturbative and background- 
independent quantization of General Relativity (GR) [1]. Interestingly, it has now been 
demonstrated that different approaches, based on canonical quantization of GR, on 
covariant quantization of GR and on formal quantization of geometry lead to the very 
same LQG theory. Although this is rather convincing, a direct experimental probe is 
still missing. One can easily argue that cosmology is the most promising approach to 
search for observational features of LQG or, more specifically, to its symmetry-reduces 
version. Loop Quantum Cosmology (LQG) [2]. 

Many efforts have been devoted to the search of possible footprints of LQG in 
cosmological tensor modes (see [3]). At the theoretical level, the situation is easier in 
this case as the algebra of constraints is automatically anomaly-free. But, as far as 
observations are concerned, scalar modes are far more important. They have already 
been observed in great details by WMAP [1] and are currently even better observed 
by the Planck mission. The question of a possible modification of the primordial scalar 
power spectrum (and of the corresponding TT Ci spectrum) in LQG is therefore essential 
in this framework. 

Gravity is described by a set of constraints. However, for the (effective) theory to 
be consistent, it is mandatory that the evolution generated by the constraints remains 
compatible with the constraints themselves. This is always true if their mutual Poisson 
brackets vanish when evaluated in fields fulfilling the constraints, i.e. if they form 
a first class algebra. This means that the evolution and the gauge transformations 
are associated with vector fields that are tangent to the manifold of null constraints. 
This obviously holds at the classical level. However, when quantum modifications are 
added, the anomaly freedom is not anymore automatically ensured. Possible quantum 
corrections must be restricted to those which close the algebra. This means that, 
for consistency reasons, the Poisson brackets between any two constraints must be 
proportional to one constraint of the algebra. This article is devoted to the search 
for such an algebra for scalar perturbations. 

Our approach will follow the one developed by Bojowald et al. in [5]. There are 
two main quantum corrections expected from LQG: inverse volume terms, basically 
arising for inverse powers of the densitized triad, which when quantized become an 
operator with zero in its discrete spectrum thus lacking a direct inverse, and holonomy 
corrections coming from the fact that loop quantization is based on holonomies, rather 
than direct connection components. In [5] the authors focused exclusively on inverse 
volume corrections. Here, we extend with work to the holonomy corrections. Scalar 
perturbations with holonomy corrections have been studied in [6] . However, the issue of 
anomaly freedom was not really addressed. Recently, a new possible way of introducing 
holonomy corrections to the scalar perturbations was proposed in [7]. Although it was 
interestingly shown that the formulation is anomaly-free, the approach is based on the 
choice of the longitudinal gauge and the extension of the method to the gauge- invariant 
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case is not straightforward. In contrast, the approach developed in our paper does not 
rely on any particular choice of gauge and the gauge-invariant cosmological perturbations 
are easily constructed. 

The theory of anomaly-free scalar perturbations developed in this paper is obtained 
on a fiat FRW background, such that the line element is given by: 

ds'^ = [-(1 + 2(\))dif + ^daBdridx" + ((1 - 2i))bah + 2dadi,E)dx''dx^'\ , (1) 

where (pjijj, E and B are scalar perturbation functions. The matter content is assumed to 
be a scalar field. This will allow us to investigate the generation of scalar perturbations 
during the phase of cosmic inflation while taking into account the quantum gravity 
effects. 

Our analysis of the scalar perturbations is performed in the Hamiltonian framework 
developed in [5l|8]. As it was shown there, the background variables are {k,p, (p, vr), while 
the perturbed variables are (5/^*, (JE'f , 5^9, ^vr). The Poisson bracket for the system can 
be decomposed as follows: 

{■) ■} = {■) 'Ifc,? + {■) ■^&K,&E + {■, ■}<^,7f + {■, ■^8ip,5-K (2) 



where 
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Here, Vq is the volume of the fiducial cell and n = 8ttG. 

The holonomy corrections are introduced by the replacement k — )■ 
classical Hamiltonian. We follow the notation introduced in |9], where 



K[n] :-- 



sin(n^7A;) 
n^7 



for n e Z/{0}, 



(3) 
(4) 
(5) 
(6) 

[n] in the 
(7) 



for 



n 



0, 



for the correction function. In cases where k appears quadratically, the integer n is fixed 
to two (See [9]). In the other cases, the integers remain to be fixed from the requirement 
of anomaly freedom. The coefficient 7 is the Barbero-Immirzi parameter and p (x p^ 
where — 1/2</3<0. In what follows, the relation 



op 



[n] = [A; cos(n/i7/c) — K[n] 



(8) 



will be useful. 

The organization of the paper is the following. In Sec. [21 the holonomy-corrected 
gravitational Hamiltonian constraint is defined. We calculate the Poisson bracket of the 
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Hamiltonian constraint with itself and with the gravitational diffeomorphism constraint. 
In Sec. [3], scalar matter is introduced. The Poisson brackets between the total 
constraints for the system under consideration are calculated. In Sec. HI the conditions 
for anomaly freedom are solved and the expressions for the counter-terms are derived. 
Based on this, in Sec. [5l equations of motion for the scalar perturbations are derived. 
The system of equations is then investigated in the case of the longitudinal gauge. 
Finally, gauge- invariant variables are found and the equations for the corresponding 
Mukhanov variables are derived. In Sec. Ej we summarize our results and draw out 
some conclusions. 



2. Scalar perturbations with holonomy corrections 



The holonomy-modified Hamiltonian constraint can be written as: 



(0) ^(2)^ 
G ' '^G I 



(1) 



(9) 
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The standard holonomy corrections are parametrized by two integers si and S2- The a, 
are counter-terms, which are introduced to remove anomalies. Those factors are defined 
so that they vanish in the classical limit (/i — 0). The counter-terms could be, in 
general, functions of all the canonical variables. We however assume here that they are 
functions of the gravitational background variables only. 

In our approach, the diffeomorphism constraint holds the classical form 
1 



Dg[N-] 



K 



(10) 



In general, the diffeomorphism constraint could also be holonomy corrected, this 
possibility was studied, e.g., in [H]. However, in LQG the diffeomorphism constraint 
is satisfied at the classical level. Therefore, if LQG is to be considered as a specific 
model of LQG, the diffeomorphism constraint should naturally hold its classical form. 
Because of this, in this paper, the diffeomorphism constraint is not modified by the 
holonomies. It is worth stressing, that the classicality of the diffeomorphism constraint 
is also imposed by the requirement of anomaly cancelation. Namely, if one replaces 
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k — )■ K[?T,] in (ITOll . the condition n = would anyway be required by the introduction of 
scalar matter. In fact, the same condition was obtained for vector modes with holonomy 
corrections [9]. 

Let us now calculate the possible Poisson brackets for the constraints H^iN] and 
Dam. 



2.1. The [h'^^Dg] 



bracket 

Using the definition of the Poisson bracket ([2]), we derive: 

'h'^[NIDg[N-]\ = -Hg [SN'^dJN] + B Dg[N-] 



S 



^ ^ d'x5N\da5E1)A^ + -^ j d'x{da5N''){5E\5l)AA, (11) 



where 



and 



B = ^ [-2K[2] + ^-(1 + ag) + K[s2] + ^e] , (12) 



A = 2fc(K[5i] + ai) + a2-2K[l]2, (13) 

-4.2 = "5 - "4, (14) 

^3 = -K[l]2-p|,K[l]2-la7 

+^(-2K[2] + ^-(1 + ag) + 2K[s2] + 2^6), (15) 

^4 = ag — tty. (16) 

The functions ^i, . . . , ^4 are the first anomalies coming from the effective nature of 
the Hamiltonian constraint. Later, we will set them to zero so as to fulfill the requirement 
of anomaly freedom. This will lead to constraints on the form of the counter-terms. 

Beside the anomalies, the |if^,DG| bracket contains the —HQ[5N'^da5N\ term, 
which is expected classically. There is also an additional contribution from the 
diffeomorphism constraint B Dg[N°']. This term is absent in the classical theory. This 
is however consistent as, for /2 — >■ 0, the B function tends to zero. 

2.2. The |ijg,i7g| bracket 



The next bracket is: 

'i/g[iVi], i/g[iV2]} = (1 + a3)(l + a,)DG 

N r 

— / d^xd''{5N2-SNi){d^SKl){l + a3)A 
Jt. 



—d''(5N2-SNi) 
.P 



N 



where 



<fx{5N2 - 6Ni){d'daSE^)A6 
N f 

+— / d'x{5N2-5N^){5'^5Ki;)Aj 
N f 

+— d^x{6N2-SN,){SlSEf)As, (17) 

A = "5 - "4, (18) 

A = (1 + ag)(K[si] + ai) - (1 + «3)(K[S2] + ae) + K[2](l + as) 



A = 4K[2]|)^(K[si] + ai) - {nif + 2p^K[l]2^ M^^'"'^ + 
3 1 \ 

1 + 2«5 - j (K[l]2 + as) - 2(K[s2] + a6)(K[si] + ai) 
+2K[2](K[si] + ai), (20) 
A = ^(Kh] + «6)(K[1]' + «2) - (K[si] + «i)(K[l]2 + «7) 

+^(K[si] + ai)(K[l]2 + as) - ^K[2](K[1]2 + a^) 

+K[2]p^(K[l]2 + a2)-\ (k[1]2 + 2p|k[1]2) £(K[1]2 + c^)- (21) 

The A;---;-4.8 are the next four anomahes. Moreover, the diffeomorphism 
constraint is multiphed by the factor (1 + a3)(l + a^). 

2.3. The {Dg,Dg} bracket 

The Poisson bracket between the diffeomorphism constraints is: 

{Dg[N^],Dg[N^]} = 0. (22) 

3. Scalar matter 

In this section, we introduce scalar matter. The scalar matter diffeomorphism constraint 
is 

DMm = jjN^n{daS^). (23) 

The scalar matter Hamiltonian can be expressed as: 

[N] = Hm [N] + Hm [SN] , (24) 

where 

Hm[N] = I d'xN [(Hf + + (Hf + + , (25) 

Hm[SN] = ^ d'SN + ?^«] . (26) 



The factors in equations (1231) and (12^ are 
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(27) 



(28) 
(29) 



Here, we have introduced the counter-term aio in the factor T-i^j . Thanks to this, the 
Poisson bracket between two matter Hamiltonians takes the following form: 



aio)-D 



M 



—d''{6N2-SNi 
P 



(30) 



As will be explained later, the appearance of the front-factor (1 + aio) will allow us 
to close the algebra of total constraints. In principle, other prefactors could have been 
expected, however they do not help removing anomalies. 



3.1. Total constraints 

The total Hamiltonian and diffeomorphism constraints are the following: 

=Hg[N]+Hf,[N], (31) 
Dtot [iVT = Dg [iVT + Dm [iVi . (32) 
The Poisson bracket between two total diffeomorphism constraints is vanishing: 

{Dtot[N'^],Dtot[N^]} = 0. (33) 

The bracket between the total Hamiltonian and diffeomorphism constraints can be 
decomposed as follows: 

{Htot [N] , Dtot [iVi } = I H'^, [N] , Dtot [iVi ] + {Hg[N], Dg [iV^ ~ 

Hg[N],DM[N^]}. (34) 
The first bracket in the sum ( IMl) is given by 

H'^[N], Aot[iVli = -H^,m-dJN]. (35) 
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The second contribution to Eq. 
vanishing: 

'/7g[iV],A,[iVl} = 



is given by fllll] . while the last contributions is 



(36) 



The Poisson bracket between the two total Hamiltonian constraints can be decomposed 
in the following way: 



{Htot[Ni],Htot[N2]} 



[iVi] , [N2] I + {Hm [Ni] , Hm [N2] } 
HSm,HM[N2]} - {N, ^ N2) 



+ 



The contribution from the last brackets can be expressed as 



H2m,HM[N2] 
1 



(iVi ^ N2) 



2p3 

p2 



+ / d-'xN{5N2-6N^){6'j5K;;) 



+3 / d-'xN{5N2-SNi) 



,{(f)6ip ) Aw 



n3 



2p 



vi^)] Mil 



vr 



+- I d'xN{SN2 - SN.MSE^) ( ^ ) A2 
+^ l_d'xN{6N2-6N,){6i6E^)V{^)A,s, 



where 
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803 
dk 



Aio = K[2] - - «i, 
^11 = 



d 3 1 

^(K[si] + ai) + -(1 + as) - + «4), 



A12 = -~(K[1]2 + a2) + 5(K[si] + ai) - 5K[2] + K[s2] + ae, 



As = ^^(^[1]' + «2) + K[SiJ + «i 
The functions A, • • • , -^13 are the last five anomalies. 



K[s2] - ae- 



(37) 



(38) 

(39) 
(40) 
(41) 

(42) 

(43) 



4. Anomaly freedom 

The requirement of anomaly freedom is equivalent to the conditions Ai = for 
i = 1,...,13. 

Let us start form the condition Ag = 0. Since 03 cannot be a constant, this 
condition implies a-s = 0. The condition = gives ai = K[2] — K[si]. Using this, 
the condition Ai = 0, can be written as 02 = 2K[1]^ — 2A;K[2]. The conditions A2 = 
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and ^5 = are equivalent and lead to = a^. Based on this, the requirement = 0, 
leads to: 

l + a^ = = cos(2/i7^) =: fi. (44) 

k 

For the sake of simplicity we have defined here the f2-function. With use of this, the 
condition ^6 = leads to 

ae = K[2]{2 + ag) - K[s2] - kn. (45) 

So, equation (H2l) simplifies to 

Ai2 = agK[2]. (46) 

Therefore, requiring Au = is equivalent to the condition ag = 0. Furthermore, A4 = 
gives a-i = ag. The expression for can be derived from the condition ^3 = 0. Namely, 
using Eq. ( l56l) one obtains: 

aj = 2(2/3 - 1)K[1]2 + 4(1 - l3)kK[2] - 2Pn. (47) 

The condition = is fulfilled by using the expressions derived for ai, a2 and a^. 
The last two anomalies ( l20l) and ( |2T]) can be simplified to: 

A = 2(l + 2/3)(fiK[l]2-K[2]2), (48) 
A = ^(l + 2/3)(K[2]2-l]K[l]2). (49) 

The anomaly freedom conditions for those last terms. Ay = and = 0, are fulfilled 
if and only if /3 = -1/2. 

It is also worth noticing that the function B given by Eq. ( fT2l) is equal to zero 
when the expression obtained for is used. There is finally no contribution from the 
diffeomorphism constraint in the lif^,/}^! bracket. 

Using the anomaly freedom conditions given above, the bracket between the total 
Hamiltonian constraints simplifies to 

'N 



tot 



■d''(5N2-6Ni] 
P 



+ (aio - a4)-D 



M 



—d''(SN2 - SNi) 
P 



(50) 



The closure of the algebra of total constraints implies the last condition aio = 04 = Q—1. 



To summarize, the counter-terms allowing the algebra to be anomaly-free are 
uniquely determined, and are given by: 



ai =K[2] -K[si], (51) 

aa = 2K[1]2 - 2kK[2], (52) 

«3 =0, (53) 

a4 =n-l, (54) 

as =n-l, (55) 
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ae = 2K[2] - K[s2] - kVt, 
aj = - AK[lf + 6^K[2] - 
as = - 4K[1]2 + 6A;K[2] - 

"9 = 0, 

aio = — 1. 



(56) 
(57) 
(58) 
(59) 
(60) 

a 10 are vanishing in the 



It is straightforward to check that the counter-terms «i, 
classical limit (/i — 0), as expected. 

Those counter-terms are defined up to the two integers Si and S2, which appear in 
( ISTl) and ( l56l) . However, in the Hamiltonian (jO]), the factor ai appears with ]K[si] 
and the factor appears with ]K[s2]. Namely, we have K[si] + ai = K[2] and 
K[s2] + ae = 2]K[2] — kQ. Therefore, the final Hamiltonian will not depend on the 
parameters Si and S2- No ambiguity remains to be fixed. 

Moreover, the anomaly cancellation requires 

(3 = (61) 

which fixes the functional form of the /i factor. The fact that anomaly freedom requires 
/3 = —1/2 is a quite surprising result. The exact value of /3 is highly debated in LQC 
The only a priori obvious statement is that /3 G [—1/2, 0]. The choice (3 = —1/2 is called 
the /2— scheme (new quantization scheme) and is preferred by some authors for physical 
reasons ^U\. Our result seems to show that the /i— scheme is embedded in the structure 
of the theory and this gives a new motivation for this particular choice of quantization 
scheme. The quantity fi'^p can be interpreted as the physical area of an elementary 
loop along which the holonomy is calculated. Because, in the /i— scheme, /i^ oc p~^, the 
physical area of the loop remains constant. This elementary area is usually set to be 
the area gap A derived in LQG. Therefore, in the /2— scheme. 



/i 



(62) 



4.I. Algebra of constraints 

Taking into account the previous conditions of anomaly-freedom, the non-vanishing 
Poisson brackets for the gravity sector are: 

= - H2[5N''dJN], 



H^[N],DGm 
Hg[N,],Hg[N,[ 



G 



—d^{SN2 - 6N^] 
P 



(63) 
(64) 



This clearly shows that the gravity sector is anomaly free. The remaining non- vanishing 
brackets are: 



{HM[N],Dtot[N^]} = -HMlSN^dJN] 
{HM[m],HM[N2]} =nDM 



—d\5N2 - SNi] 
P 



(65) 
(66) 
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The algebra of total constraints therefore takes the following form: 



{DUNtlDt,t[N^]} 
{Ht,t\NlDt,t[N'']} 

{Htot[NilHtot[N2]] 



0, 

Dtot 



N 

n—d''(6N2 
. p 



(67) 
(68) 

(69) 



Although the algebra is closed, there are however modifications with respect to the 
classical case, due to presence of the factor Q in Eq. ( 169|) . Therefore, not only the 
dynamics, as a result of the modification of the Hamiltonian constraint, is modified but 
also the very structure of the space-time itself is deformed. This is embedded in the 
form of the algebra of constraints. The hypersurface deformation algebra generated by 
( l69l) is pictorially represented in Fig. 14.11 As Q G [—1,1], the shift vector 

N 

iV« = -cos{2fi-fk)d''(6N2 - 5Ni) 
P 




Figure 1. Pictorial representation of the hypersurface deformation algebra (|69p . 



= n—d''(SN2 - SNi) (70) 
P 

appearing in fl69l) can change sign in time. 

In order to see when this might happen let us express the parameter Q as: 

n = cos(2fi-fk) = 1-2—, (71) 

Pc 

where p is the energy density of the matter field and 

= ^ = A- (72) 

In the low energy limit, p — 0, the classical case {Q — ?■ 1) is correctly recovered. 
However, while approaching the high energy domain the situation drastically changes. 
Namely, for p = Pc/2, the shift vector flTUl) becomes null. At this point, the maximum 
value of the Hubble parameter is also reached. The maximum allowed energy density 
is p = Pc and corresponds to the bounce. Then the shift vector fl7UI) fully reverses with 
respect to the low energy limit. One can interpret this peculiar behavior as a geometry 
change. Namely, when the universe is in its quantum stage (p > Pc/2), the effective 
algebra of constraints shows that the space is Euclidian. At the particular value p = y, 
the geometry switches to the Minkowski one [11]. This will become even clearer when 
analyzing the Mukhanov equation in Sec. O The consequences of this have not yet 
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been fully understood, but it is interesting to notice that this model naturally exhibits 
properties related to the Hartle-Hawking no-boundary proposal |12j . 

5. Equations of motion 

Once the anomaly-free theory of scalar perturbations with holonomy corrections is 
constructed, the equations of motion for the canonical variables can be derived. This 
can be achieved through the Hamilton equation 

f = {f,H[N,N'^]}, (73) 

where the Hamiltonian H[N, N""] is the sum of all constraints 

H[N, N-] = h2[N] + Hm[N] + Dg[N-] + Dm[N% (74) 

5.1. Background equations 

Based on the Hamilton equation (173!) . the equations for the canonical background 
variables are the following: 



(75) 



2p3 

p = 2iVVpK[2], (76) 

^ = (77) 

^ = -Np^/^V^{^). (78) 

In the following, we choose the time to be conformal by setting N = y/p. The " ■ " then 
means differentiation with respect to conformal time rj. 

Eqs. ( 1771) and ( 1781) can be now combined into the Klein-Gordon equation 

ip + 2K[2]^ + pV^{^) = 0. (79) 

Eq. (1761) . together with the background part of the Hamiltonian constraint 
1 dH 1 



v,3N i-^-'^^mm+p 

lead to the modified Friedmann equation 



3/2 



7r2 



2f 



0, (80) 



3 V 

Another useful expression is: 



n' = p^p[i-^]. (81) 



7r2 



m[lf = — +pV{^). (82) 
Here % stands for the conformal Hubble factor 



The energy density and pressure of the scalar field are given by: 



TT 



For the purpose of further considerations, we also derive the relation 

= kK[2] - it, 



2p2 

which comes from Eq. fl75|) combined with f lHOj) . 
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(86) 



5.2. Equations for the perturbed variables 
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The equations for the perturbed parts of the canonical variables are: 
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Vp 
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5.5. Longitudinal gauge 

As an example of application we will now derive the equations in the longitudinal gauge. 
In this case, the E and B perturbations are set to zero. The line element ([1]) therefore 
simplifies to 

ds^ = [-{1 + 20)c/r/2 + (1 _ 2i/j)6abdx''dx''] , (91) 



14 



where (f) and are two remaining perturbation functions and a is the scale factor. From 
the metric above, one can derive the laps function, the shift vector and the spatial 
metric: 

av/1 + 20, (92) 
0, (93) 
2i^)Sab. (94) 



N 

Qab 



The lapse function can be expanded for the background and perturbation part as 
N = N + 6N, where 

N =VP = a, (95) 
SN = N(f). (96) 

Using Eq. ( 194|) . the perturbation of the densitized triad is expressed as: 

6Et = -2p^St. (97) 

The time derivative of this expression will also be useful and can be written as: 

6E^ = -2p{2K[2]ij + ^)51. (98) 

Let us now find the expression for the perturbation of the extrinsic curvature SK^^ 
in terms of the metric perturbations (f) and ip. For this purpose, one can apply the 
expression Wl\ to the left hand side of (1571) . The resulting equation can be solved for 
5Kl^, leading to: 



5Ki = -5i^(^^ + knij + W(2]<p). 
The time derivative of this variable is given by 



(99) 



5Ki = 61 



- (/.K[2] - K[2](j) 



ip — kVtip + if) 



n 



kn + 0K[2] 



n 



(100) 



Applying fllOOp to the left hand side of flHHjl . the equation containing the diagonal part 
as well as the off-diagonal contribution is easily obtained. The off-diagonal part leads 
to 

dad\<j)~tlj) = 0. (101) 

This translates into ip = (p. In what follows, we will therefore consider only. The 
diagonal part of the discussed equation can be expressed as: 



+ 



+ 



+ 2K[2]2-K[2]|^ 



One case now use the diffeomorphism constraint 



(102) 



K- 



pd,{6tdK'^)-p{dkSK'^)-kS'^id,SE^,) + K7T{dM = 0.(103) 



15 



With the expressions for SK^^ and SEf, it can be derived that 



0K[2] = AirGn^dcSif. 



(104) 



The next equation comes from the perturbed part of the Hamiltonian constraint: 



6H[N, N"] 
6{6N) 



1 

2k 
2 



-4./pK[2]6'j6Ki - — (3K[1]2 - 2kK[2]) 6i6E. 



VP , 



non 



IT 



P 



;3/2 



2p3/2 

6i6ET 



2p 



2p 



0. 



(105) 



Using the expressions for 6K^^ and SE^, this can be rewritten as: 



K[2] + 2K[2]^ 



4nGn [^6<^ + p6<^V^{^)] . (106) 



fiV^(/) - 3K[2]0 
The last equahty comes from fl89l) and fl90l) : 

6<p + 2K[2]5v3 - nV^S^ + py^^((^)(5(^ + 2pV^^((^)(/) - 4(^0 
The equations fll02p . (11041) and (I106P can be now combined into: 



0. 



+ 2 



n 



with the quantum correction 



+ e 



P 



(107) 



(108) 



(109) 



and the squared velocity 



(110) 

The squared velocity of the perturbation field </> is equal to f2. Because — 1 < < 1, 
the speed of perturbations is never super-luminal. However, for f2 < perturbations 
become unstable (c^ < 0). This corresponds to the energy density regime p > y, where 
the phase of super-inflation is expected. 

At the point p = y, the velocity of the perturbation field (f) is vanishing. Therefore, 
perturbations don't propagate anymore when approaching p = y, where the Hubble 
factor reaches its maximal value. Moreover, at this point, the quantum correction 
e — > oo. Because of this, Eq. f llOSp diverges and cannot be used to determine the 
propagation of the perturbations. However, as shown in the next section, the equation 
for the gauge- invariant Mukhanov variable does not exhibit such a pathology. 

It is interesting to notice that the equations of motion derived in this subsection 
are the same as those found in [7j. This is quite surprising, because they were derived 
following independent paths. In our approach, we have introduced the most general 
form for the holonomy corrections to the Hamiltonian. Then, by adding counter-terms, 
anomalies in the algebra of constraints were removed. On the other hand, the method 
proposed in [7] is based on the diagonal form of the metric in the longitudinal gauge. 
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This enables one to introduce holonomy corrections in almost the same way as in the 
case of a homogeneous model. It was then shown that a system defined in this way stays 
on-shell, that is, is free of anomalies. The non-trivial equivalence of both approaches 
may suggest uniqueness in defining a theory of scalar perturbations with holonomy 
corrections in an anomaly-free manner. 

5.4- Gauge invariant variables and Mukhanov equation 

Considering the scalar perturbations, there is only one physical degree of freedom. As it 
was shown in [13], this physical variable combines both the perturbation of the metric 
and the perturbation of matter. The classical expression on this gauge- invariant quantity 
is: 

V = a{r]) ( 6^^' + , (111) 



and its equation of motion is given by 

V -V^v --V = 0, (112) 

z 

where 

z = aiv)^. (113) 

In the canonical formalism with scalar perturbations, the gauge transformation of 
a variable X under a small coordinate transformation 

x^^x^ + e^ ; = (^o,9'^0, (114) 

is given by (see [8] for details): 

%o,^]X={x,i/(2)[7ve°] + (115) 

and it is straightforward to see that, classically. 

This means that v is diffeomorphism-invariant and can be taken as an observable. 

Taking into account the holonomy corrections introduced in this paper, the Q 
function will modify the gauge transformations of the time derivative of a variable X, 
so that 

- = ^ • ^[0,«0]^- (117) 

Using this relation and gauge transformations of the metric perturbations 

= - me, (118) 

%o,5]0 =e° + K[2]e°, (119) 
= e, (120) 

= i (121) 
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one can define the gauge-invariant variables (Bardeen potentials) as: 

^ =ij-^{B-E), (123) 

6^^' = 6^ + ^{B~E). (124) 

The normalization of these variables was set such that, in the longitudinal gauge 
{B = = E), we have $ = 0, \1/ = t/^ and 6ip'^^ = 6ip. It is possible to define the 
analogous of the Mukhanov variable fill II) : 

v:=VP i^V^' + ) • (125) 



c^V^^ = (126) 



5(^^^+2K[2](5v9^'^-fiV2<5<^^^+py^^((^)<5(p^^+2p\/^((^)^-4(^^^^ = 0, (127) 



K[2] 

Writing the equations for \& and 6ip^^ , which are 



and 



^ + 2 






^ + 2 



















one obtains equation for the variable (11251) : 

V - nV'^v --v = 0, (128) 

z 

which corresponds to the Mukhanov equation for our model. As we see, the difference 
between the classical and the holonomy-corrected case is the factor fl in front of the 
Laplacian. This quantum contribution leads to a variation of the propagation velocity 
of the perturbation v. This is similar to the case of the perturbation considered in 
the previous subsection. The main difference is that there is no divergence for p = Pc/2 
and the evolution of perturbations can be investigated in the regime of high energy 
densities. It is once again worth noticing that for p > Pcf^, ^ becomes negative and 
Eq. (I128p changes from an hyperbolic form to an elliptic one. This basically means that 
the time part becomes indistinguishable from the spatial one. This can be interpreted 
as a transition from a Minkowskian geometry to and Euclidean geometry, as mentioned 
earlier. 

Finally, it is also possible to define the perturbation of curvature TZ such that 

7^ = ^. (130) 

Based on this, one can now calculate the power spectrum of scalar perturbations. This 
opens new possible ways to study quantum gravity effects in the very early universe. 
Promising applications of the derived equations will be investigated elsewhere. 
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6. Summary and conclusions 

In this paper, we liave investigated the theory of scalar perturbations with holonomy 
corrections. Such corrections are expected because of quantum gravity effects predicted 
by LQG. They basically come from the regularization of the curvature of the connection 
at the Planck scale. Because of this, the holonomy corrections become dominant in the 
high curvature regime. The introduction of "generic type" holonomy corrections leads 
to an anomalous algebra of constraints. The conditions of anomaly freedom impose 
some restrictions on the form of the holonomy corrections. However, we have shown 
that the holonomy corrections, in the standard form, cannot fully satisfy the conditions 
of anomaly freedom. In order to solve this difficulty, additional counter-terms were 
introduced. Such counter-terms tend to zero in the classical limit, but play the role 
of regularizators of anomalies in the quantum (high curvature) regime. The method 
of counter-terms was earlier successfully applied to cosmological perturbations with 
inverse-triad corrections [5]. 

We have shown that, thanks to the counter-terms, the theory of cosmological 
perturbations with holonomy corrections can be formulated in an anomaly-free way. 
The anomaly freedom was shown to be fulfilled not only for the gravity sector but also 
when taking into account scalar matter. The requirements of anomaly freedom were 
used to determine the form of the counter-terms. Furthermore, conditions for obtaining 
and anomaly-free algebra of constraints were shown to be fulfilled only for a particular 
choice of the /2 function, namely for the /i— scheme (new quantization scheme). This 
quantization scheme was shown earlier to be favored because of the consistency of the 
background dynamics [10]. Our result supports these earlier claims. 

In our formulation, the diffeomorphism constraint holds its classical form, in 
agreement with the LQG expectations. The obtained anomaly-free gravitational 
Hamiltonian contains seven holonomy modifications. It was also necessary to introduce 
one counter-term into the matter Hamiltonian in order to ensure the closure of the 
algebra of total constraints. There is no ambiguity in defining the holonomy corrections 
after imposing the anomaly-free conditions. The only remaining free parameter of the 
theory is the area gap A used in defining the fl function. This quantity can however be 
possibly fixed with the spectrum of the area operator in LQG. Based on the equations 
derived in this paper it will also be possible to put observational constraints on the value 
of A and, hence, on the critical energy density pc- 

Based on the studied anomaly-free formulation, equations of motion were derived. 
As an example of application, we studied the equations in the longitudinal gauge. We 
have also found the gauge-invariant variables, which are holonomy-corrected versions 
of the Bardeen potentials. Using this, we have derived the equation for the Mukhanov 
variable. This equation can be directly used to compute the power spectrum of scalar 
perturbations with quantum gravitational holonomy corrections. Similar considerations 
were studied in the case of inverse-triad corrections jT4]. In that case, observational 
consequences have been derived and compared with CMB data [15l |16] . 
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